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Expectation values of the electromagnetic field and the electric current are introduced at space- 
time resolution which belongs to the quantum domain. These allow us to approach some key features 
of classical electrodynamics from the underlying QED. One is the emergence of the radiation field in 
the retarded solution of the Maxwell equation, derived from an action principle. Another question 
discussed is the systematic derivation of the polarizability of a charge system. Furthermore, the 
decoherence and the consistency of the photon field are established by a perturbative calculation of 
the reduced density matrix for the electromagnetic field within the Closed Time Path formalism. 
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I. INTRODUCTION 

Classical electrodynamics consists of a set of equations of motion for the electromagnetic (EM) field and the electric 
current, the Maxwell and the mechanical equations, respectively. But the classical theory should be derived from 
QED rather than used as a starting point to quantize the EM field. When approached in this manner, these equations 
become approximate, govern the expectation values of the EM field and the electric current, and are the subject of 
an infinite hierarchy of quantum corrections. 

The degrees of freedom common in quantum and classical physics are the expectation values of local fields. But we 
have to keep in mind that such expectation values are not necessarily classical. The crucial quantity which controls 
the classicality of the field expectation value is the number of degrees of freedom N = n£ 3 within the internal distance 
scale of the field, the correlation length £ of a many-body system at particle density n. The UV cutoff A, the maximal 
energy scale of QED, will be kept large but finite in this work: therefore, the internal scale of the field expectation 
value £ > 1/A may fall into the quantum or classical domain, depending on initial conditions and external sources. For 
N ~ 1 the local expectation value does not follow the laws of classical physics and will be called a subclassical field. 
Such fields are used for instance in the density functional method [5] . These fields are useful to monitor the quantum- 
classical transition because they bridge the gap between observables and expectation values, used on the microscopic 
and macroscopic domain, respectively. This can be seen simplest within the framework of the renormalization group 
where the classical action governing the dynamics of the subclassical field can be set as an initial condition for the 
integration of the renormalization group equations in the macroscopic domain. This point of view motivated this work, 
the study of some necessary ingredients of the quantum-classical transition by means of the dynamics of subclassical 
fields. It includes the study of the dynamical origin of decoherence [1, 2] and consistency [3], the breakdown of the 
time reversal invariance [4] and the way the expectation values are formed at the microscopic scale. 

A reassuring feature of subclassical fields is that their dynamics is stable and well defined in the UV. In fact, 
the particles can not be localized with a precision better than their Compton wavelength therefore the subclassical 
fields are smeared out at shorter distances. This scale is independent of the underlying quantum field theory. The 
subclassical field is classical in the sense that quantum fluctuations are averaged out. But it is still in the quantum 
domain due to its fine spatial resolution because the decoherence within this space region is not strong enough to 
restore classical probabilities, in particular, the additivity for exclusive events. 

The present work is the continuation of the program outlined in Rcf. [20] where the equations of motion for 
the expectation value of the EM field and the electric current were discussed in QED. We go beyond that work 
and render that approach more realistic by including a reservoir of charges at finite density and by calculating 
explicitly and by interpreting the vacuum polarization effects in the equation of motion for the EM field and in 
decoherence/consistency. The loop expansion will be employed to include the vacuum polarization effects and we 
shall be satisfied by calculating the leading order, one-loop approximation to the linearized variational equations of 
motion. The current terminology of the traditional construction of the averaged fields in classical electrodynamics can 
be borrowed from the renormalization group method, it is a special blocking. In the present work we are satisfied by a 
simpler realization of this idea, the coarse graining generated by the elimination of the charged degrees of freedom from 
the dynamics without lowering the space resolution. The decoherence and consistency become simple perturbative 
effects in this setup, they are generated by the collective modes of the many-body vacuum. The dynamical breakdown 
of the time reversal invariance can be traced back to the collective modes as well, it is driven by the fmiteness of 
the life-time of the quasiparticles. Both the decoherence/consistency and dynamical breakdown of the time reversal 
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invariance are actually generated by the same piece of the photon polarization tensor. This is how the dynamical 
breakdown of the time reversal invariance and decoherence appear coupled in the classical limit. Note that there is 
no need of a particular definition of pointer states in this procedure, quasiparticles arising form our coarse graining 
provide a natural, dynamically generated robust basis in the Fock space for the study of the reduced density matrix. 

The breakdown of the time reversal invariance brings the boundary conditions in time into the foreground. The 
other goal of this work is the establishment of the dynamical symmetry breaking and its clear separation from the 
explicit breakdown by the boundary conditions in time. The not entirely trivial issue of the latter is the way the 
external sources rearrange the radiation field in the variation principle to comply with the given boundary conditions. 
The simplest way to regulate the corresponding formal, 0/0 terms in the equations of motion is the introduction of 
the usual ie term in the quadratic part of the action in classical field theory, as well as in its quantum counterpart. 

The quantum corrections to the space-time dependence of the averages of local observables are usually calculated 
by means of the linear response formalism. We need a more powerful method, capable of dealing with the eventual 
nonpcrturbative issues of the quantum-classical crossover regime to arrive at a reliable description of quantum effects, 
turning classical electrodynamics into QED at short distances. The Closed Time Path (CTP) method developed by J. 
Schwinger [6] provides us with a nonperturbative setting not only for the determination of the space-time dependence 
of local observables but for other key ingredients of the quantum-classical transition, to the reduced density matrix 
and the possible dynamical breakdown of the time reversal invariance, too. 

The CTP method has been extended a number of times [7-15]. Despite being the only way of obtaining expectation 
values as opposed to transition amplitudes reproduced by conventional quantum field theory its inherent algebraic 
complications somehow limited its applications. The generator functional for connected Green functions, our starting 
point given by Eq. (7) below is closely related to the influence functional introduced in Ref. [16]. One may simple 
say that all that is presented below is the application of the influence functional scheme to derive the effective action 
and the equations of motion for the expectation value: of local operators. The: solution of the latter provides us with a 
partial resummation of the perturbation series with rich physical content. Furthermore, certain details of the influence 
functional contain decoherence/consistency and irreversibility, the key conditions for the classical limit. The reduced 
density matrix and dissipative forces have been studied in the framework of nonrelativistic quantum mechanics in a 
manner reminiscent of the CTP formalism [17]. The unusual feature of quantum transition amplitudes, the possibility 
of fixing the initial and the final states of the motion has already been addressed in a manner, equivalent to the 
CTP formalism in Ref. [18]. The Schwinger-Dyson equation, the removal of the UV divergences, the equation of 
motion for the EM field derived in the CTP formalism and the pair creation generated by a homogeneous external 
electric field have already been thoroughly studied for a large number of charged fermions in Ref. [21]. Our results 
are complementary since they retain the adiabatic, IR tail of the vacuum polarization only. It is worth mentioning 
that the problem of establishing variational equations for the expectation values [19] can be resolved by a suitable 
parametrization of the degrees of freedom [20] . The role of collective modes in the quantum-classical transition can 
conveniently be studied within the CTP scheme [22], as well. 

We start in Section II with a short introduction of the CTP formalism, including the systematic construction of 
the Heisenberg representation in terms of Green functions. Instead of following the historical route, the motivation of 
this scheme by the calculation of the expectation values rather than transition amplitudes, an alternative argument is 
presented. The CTP formalism is motivated by the generalization of the classical action principle for initial condition 
problems which produce the retarded solutions. Beside the CTP case we need the Open Time Path (OTP) version, as 
well, to determine the reduced density matrix. Some functionals, used subsequently are introduced here, too, namely 
those which generates the connected Green functions and their Legendre transform, the effective action. 

A rather general issue is addressed in Section III where a nonrenormalization theorem is shown for the Maxwell 
equations when all charges are included. One finds dressing only when charges are ignored. This result is in agreement 
with the usual derivation of the dielectric constant and the magnetic permeability in classical electrodynamics. 

The double-time formalism is introduced for QED in Section IV, followed by the calculation of the quadratic part 
of the generator functional for the connected Green functions and the effective action for the electric current and the 
EM field, as well as a short discussion of the ways the mass-shell modes of the photon and the zero modes of the 
free EM field action satisfy the desired boundary conditions in time by formal, 0/0- type terms in the equations of 
motion. The one-loop result for this functional is summarized in Section V. 

The explicit calculation of the reduced density matrix for the EM field is carried out in Section VI. The real part 
of its logarithm controls the dielectric and the magnetic permeability functions and determines the quasiparticles, the 
imaginary part is responsible for the finite life-time of the quasiparticles, and the decoherence and the consistency of 
the EM field. Finally, Section VII summarizes our conclusions. 

The Appendixes contain some technical details. A short summary of the CTP propagators in the vacuum is 
presented in Appendix A, followed by the extension to thermodynamical equilibrium in Appendix B. 
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II. TIME ARROW AND THE DOUBLE-TIME FORMALISM 



The dynamical breakdown of the time reversal invariance is a provocative, not yet understood problem [4] . Weak 
interactions set apart, the fundamental forces and dynamical equations of Nature are invariant under a reversion of 
the direction of time. Nevertheless, daily experiences show that external perturbations influence the systems at later 
times only, and therefore the retarded solutions of the time reversal invariant equations of motion should be used. We 
set up our quadrature to solve the differential equations accordingly, by specifying the initial conditions. The resulting 
scheme obviously does not obey time reversal invariance. Such an explicit breakdown of the symmetry is a trivial 
phenomenon, the real, dynamical issue is the simplicity of the initial conditions, and the absence of the radiation field 
in electrodynamics or excited states, in general, in quantum mechanics. 

Even this trivial issue poses an interesting question. How can such an initial condition problem be dealt with 
in physics? The Newton equation, being a second order differential equation, requires two input data per degrees 
of freedom to fix a solution uniquely. For a differential equation we can choose freely among parameters of the 
motion, either at the initial or at the final time. This allows us to treat problems where the initial or the final 
condition is simple. A typical example is the use of the retarded or advanced Lienard-Wicchcrt potentials in classical 
electrodynamics. But the fundamental, most flexible layer of classical mechanics is the variational principle where 
boundary problems can be defined only because the variation performed at the end point of the trajectory in an initial 
problem cancels the generalized momentum. Such a shortcoming of the variation method is not a serious problem 
in classical dynamics since we can switch back to the canonical equations of motion anytime and set up the desired 
Cauchy problem. The situation is more involved in quantum mechanics, where the action principle remains important. 
In fact, after having set up a problem in the operator formalism, the actual calculation of the transition amplitudes 
or expectation values is usually done in the path integral formalism. Hence the goal of generalizing the variation 
principle to cope with retarded or advanced forces remains an important question. 

How can we construct a variational method which automatically produces the retarded solution of the equations 
of motion? One should modify the procedure in such a manner that it is enough to provide the coordinates and the 
velocities at t = ti, the beginning of the time evolution. But what does one do with the variation of the trajectory at 
the final time t = tfl A simple way to render the variation at the final time consistent with the classical trajectory 
is to cancel altogether its contribution to the equation of motion. This can be achieved by extending the dynamical 
problem for twice as long as < tcrp < 2T = 2(tf — ti). We follow the usual trajectory x + (t) for the time interval T 
then we flip the time arrow and follow the time backward along another trajectory x~(t), which returns to the initial 
conditions. The closed trajectory obtained in this manner is 



The variation problem for such a time reversed trajectory is well defined, and it yields a trivial equation of motion 
0=0 for t = tf and produces the retarded solution. The quantum analogy of this scheme is well known, it is the CTP 
formalism of Schwinger [6]. We shall need a slight generalization of this scheme, the OTP fromalism [23] to access 
the reduced density matrix of subsystems to cope with a genuine quantum effect, the decoherence. 

The CTP formalism was originally developed as the standard perturbation expansion for expectation values in 
the Hcisenberg representation of Quantum Field Theory [6]. When Feynman's insight to QED, based on a simpler 
scheme of transition amplitudes became widely accepted, then the perturbation expansion for the expectation values 
was gradually abandoned. This shift of attention seemed to be further justified by the simplicity and generality of the 
path integral representation of the transition amplitudes. The interest in the CTP formalism was kept in condensed 
matter physics, where expectation values corresponding to given, time-dependent external perturbation are important. 

The expectation value 



obtained in the Heisenberg representation, where where An(t) — e'(* Ae '(* ti ^ H and |^>(i»)) = can be 
written as 



in terms of the initial density matrix pi — |^>j)(^j|. This expectation value can be related in a simple manner to the 
transition amplitude 




(1) 



(2) 



{i){t)\A\^(t)) = TrlAe-^-^Ve^"*^]. 



(3) 



A = {ipile-^f-^Ae-^-^ 



(4) 
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only when the initial state is an eigenstate of the Hamiltonian. This condition is not met in a realistic situation 
where the initial state contains complicated, collective excitations. Expectation values like (3) will be obtained in 
the CTP formalism, used in this work to derive the effective action whose Euler-Lagrange equation is satisfied by the 
expectation values of the EM field and the electric current. 



A. Generator functional for connected Green functions 



We recall first some general features of the CTP formalism and next we introduce this scheme for QED. The 
expectation value (3) motivates the introduction of the generator functional 

e iW[j+,j-] = IVT^-*/*/ dt f d. 3 x[H(x)-j+(x)0(x)]^ T ^-i ft/ dt f d 3 x[H(x)+j- (x)O(x)]^ ^ 

where j^ 1 are two sources coupled to a local operator O and we use the units h = c = 1. This generator functional can 
easily be given in the path integral formalism. Indeed, the time evolution operator to the left of pi has the standard 
path integral representation. The same holds for the time evolution operator on the right except that the replacement 
j + — > — j~ is carried out, as well. Note that the Hermitian conjugation of the operator standing to the right of pi 
leads to anti-time ordering. Consequently, the number of degrees of freedom is doubled, cf> + and </>~ will be used to 
denote the CTP doublets, field trajectories for the left and right time evolution operator, respectively. We therefore 
have pairs of trajectories in the full path integral 



e «vU + J-]= f D[^]D[#]D[^r]tf*[^-]* o [0+] /"**' 0/ C [0+]^[0-] e ^[0 + ]+^cT[0+]-,5[0-]- J 5cT[0-]+ 4J +o + + y -o- 



(6) 

where is constructed by means of the field </> ± , the integral measure £ ) [0 ± ] includes the field variables for tj <t <tj 
and the space coordinate of <fi is suppressed. The initial condition pi — |0)(0| is used where |0) stands for the 
noninteracting vacuum state with wave functional \I/o[0] and the trajectories (f> + and 0~ are joined at the final time 
due to the trace operation in Eq. (5). The boundary conditions in time do not influence UV divergences therefore 
the counterterms contained in Sct[(/)} are imposed separately for each time axes. It is not difficult to check that 
the boundary conditions in time, used in Eq. (6) can be incorporated by the replacement <S'[^ + ] — S[<j)] = 

S[4> + ] — S{4>~] + Sbc[4> + i 4>~] where Sbc 1S a quadratic expression of the fields taken at t = t j and t = tf. The scalar 
product fg of space-time functions stands for space-time integration, fg = J d 4 xf x g x = J x f x g x , vector and CTP 
indices being summed, if necessary, eg. jA = X)<t=± I x Jx^^x f° r the minimal coupling. Note that the expectations 
values, calculated for unitary time evolution j + = —j~ are independent of the value of this final time tf as long as 
they are considered before tf. The bonus of these apparent complications in the CTP formalism is a way to handle 
the dynamical breakdown of the time reversal invariance, decoherence and dissipative forces [20, 22]. 

We consider in this work elementary fermions of electric charge e„, n being a flavor index, interacting with the EM 
field. For this kind of system, the quantities like (3) can be obtained by the use of the generator functional 



e iW[, 



r [a+,a-,j+,j-] = rp r T[e~ i ^i dt S d 3 x[H{x)-^ n a+{x)J n {x)-j+{x)A{x)]^ pSr *^ e i f t ' dt J d 3 x[H(x)+Y. n a~ (x)J n (x)+j~ 

where T* denotes anti-time ordering, are two sources used to generate the EM field A and a are two vectorial 
sources coupled to the EM current -ipj^ip = J^. Eq. (7) can be written as 



:^'J = J D^n^p^^EnE.^^ (8) 



using the notation 



*■-(£)■ *-(£)• *•-($■ ^e r >- 



In Eq. (8), the CTP inverse-propagators 
fi-i _ fip-m n +ie V^-i A-i _ f^T + £DL + ie \ , rifV , 
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already take into account boundary conditions in time, represented by the full CTP matrices G B q and D B c- The 
notation 

f)aab 

T ab = g ab _ L ab^ La b = ^ 

will be used for the transverse and longitudinal projection operators. The counterterm Sct will be suppressed below 
because it is given by the usual transition amplitude formalism. 

The CTP propagators are introduced in Appendix A for free fields. The bosonic two-point functions, such as the 
photon propagator, possess the structure 

" / D D+-\ (D n -Df\ (I l\ 

= D-) = { D f -D» ) +iD \ll) W 

involving three real space-time dependent functions. It is an important property of radiative corrections that the 
self energies appearing in the quadratic generator functionals display the same structure and allows us to define the 
interactive near and far field propagators. The actual calculation of the propagator for photons reveals that D n 
and D f are the near and far field Green functions of classical electrodynamics, respectively. These propagators are 
usually introduced by splitting the Lienard-Wiechert potential of classical electrodynamics into an C(r _1 ) part and 

an 0(r~ 2 ) part and the retarded and advanced propagators are defined as D a = D n ± . The propagator D is 

symmetric with respect to the exchange of its indices, cf. Eq. (Al). Therefore, the near and far field propagators 

f f 

can be identified by the symmetric and antisymmetric part of the retarded Green function, D" b — D ba , D ab = —D ba . 
Another, more microscopic separation of these two Green functions, based on the quantum level, is mentioned at the 
end of this section. Finally, we shall find in Section IV C an additional, classical way to define the near and far field 
propagators. 



B. Expectation values of the EM field and electric current 

Notice that the unitarity of the time evolution assures that expectation values can be calculated in two different 
ways, e.g. 

SW SW 

{A) = -rrf =^z (13) 

(IT |j=a=0 °J \j=a=0 

and the final time tf plays no role as long as the expectation value is considered before it, x° < tf, in the equations 
above. Such a noncompact R symmetry, called the CTP symmetry, suggests the re-parametrization [20] 

a ± = ^(l±K)±a, ^ = ±$.±$±3 (14) 
of the external sources because we have 

°J \j=a=0 

for an arbitrary real value of k. The book-keeping sources a and j must be canceled after calculating the functional 
derivatives in order to recover the physical expectation values. On the contrary, the sources a and j can be nonvanishing 
for unitary time evolution and represent the devices which are supposed to drive the system adiabatically from the 
vacuum at t = — oo to the desired initial state at t — ti. Eq. (15) shows that the fields coupled to a and j can be 
identified by the physical fields because their averages reproduce the field expectation values. We shall see later that 
the field coupled to the sources a and j will be responsible for decoherence. 

The obvious difficulty of the construction of variational equations for the expectation values is the fact that we 
have twice as many degrees of freedom in the CTP formalism as in the classical counterpart. This problem can be 
overcome by considering the effective dynamics obtained by the elimination of a degree of freedom by its equation of 
motion from each CTP doublet. But the preceding argument reveals an apparent problem in this plan. On the one 
hand, the reduplication of the degrees of freedom in the CTP formalism, A — > A ± 7 produces two identical looking 
variables as far as expectation values are concerned and it is not clear which one to eliminate. On the other hand, 
decoherence, a mechanism not directly encoded by expectation values is related to the difference A d = A + — A~ only. 
There are different dynamical issues mixed in these variables. Therefore, it appears to be appropriate to go over the 
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Kcldysh basis, (A+,A~) -t (A,A d ) = ((A+ + A~)/2, A+ - A~) and eliminate A d . But there is a new problem which 
arises in this manner. In fact, it is not difficult to see that the dynamics of the field variable A and its canonically 
conjugated momentum II = is contained in the dependence of A or A d , respectively of the density matrix 

/?L4 + , A~\. The complete elimination of A d removes all information about the momentum in the effective theory. But 
the time dependence of the expectation value (A) can be used to reconstruct the expectation value (II) = do (A) . Thus 
A d should appear in the combination of A + and A~ to be retained in the effective theory. This condition is realized 
by the restriction K^Oon the otherwise free parameter [20] . 

The effective action, a device used to derive the equation of motion for the expectation values, is usually the 
functional Legendre transform of W . This step must be slightly modified because the latter object is complex in the 
double-time formalism. Since W = for the physical case with unitary time evolution we define the effective action 
as the Legendre transform of the real part of W, 

$tW[a,j] = T[j,A] + U + j A, (16) 

and 

*-f. 

by using the external sources and j^. The inverse transformation, Eq. (16) and 

—■D' >=-ji' (18) 

appear as variational equations of motion satisfied by the fields J and A. It is advantageous to find the effective 
action for the expectation values. To this end we write 

T[J,J,A,A] =mV[a,j]-aJ - aJ '- jA-jA (19) 

by means of the parametrization (14) and 

sstw - smv „ smv 7 smv , . 

J — < J — ' A = -8T> A -—- (20) 

The inverse transformation, 

sr sr . sr , sr 

a = ~Tr fl = "w' j = j = sa> (21) 

represents the Euler-Lagrange equations. The fields J and A stand for the expectation values and the auxiliary fields 
J and A incorporate the effects of quantum fluctuations. 

The imaginary part QW, left out from the construction of the our effective action is important in establishing 
decoherence. In fact, decoherence stands for the suppression of the off diagonal elements of the density matrix in 
the pointer representation. By assuming that the field variables J and A are good pointer variables decoherence 
corresponds to the suppression of the absolute magnitude of e %w ^ a '^ as J and A are increased. Such a suppression 
comes from 3iy[o,j]. 

Finally, one can construct an effective action for the physical expectation values only, 

T[J,A]=$tW[a,a,j,j}-aJ-jA (22) 

where 

snw „ smv 



The inverse transformation, 



a = -Jj' J = -5A> (24) 



give the Euler-Lagrange equations for the expectation values only. 
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C. Reduced density matrix 

We have surveyed so far the means to find equations of motion for expectation values. The right-hand side of Eq. 
(3) is not sufficient for our other goal, for the calculation of the reduced density matrix. The reduced density matrix 
for the EM field is obtained by eliminating the charged degrees of freedom and can be written as 

P [AV>,aW] = (A^\Tv ch [Ae-^-^ H Pl e^-^ H }\A^) (25) 

where the trace operation is over the Fock space of the charged particles and the state \A) is the eigenstate of the 
photon field operator with field configurations A^x) in the functional Schrodinger representation. The generator 
functional of Eq. (8) for p[A {1 \ A^} will be given in the OTP formalism. It consists of a path integral where the 
charged fields have a closed time path and the integration over the photon field is restricted to trajectories with open 
end points A+(t f ,x) = A$\x) and A~(tf,x) = A^\x). 

The OTP formalism displays entanglement of subsystems in a specially clear manner. Let us consider for instance 
QED as a closed system, described by a factorizable, pure state at the initial time. The interactions generate 
entanglement between the charges and the EM field; the charges considered alone will be found in a mixed state and 
their reduced density matrix can be obtained by eliminating the EM field. The contributions of the perturbative 
integration over the EM field configurations are labeled by graphs. Those which contain photon lines connecting the 
two time axes represent mixed state contributions. In fact, the integration over the momentum of these lines produces 
the sum of pure state contributions to the reduced density matrix. The structure (12) of the photon propagator 
reveals that the near field is responsible for the interactions within a time axis. These elementary processes contribute 
to the self energy and preserve the purity of the charge states. The far field connects the two time axes and represents 
the mixing resulting from the entanglement of the charge-EM field subsystems. 

The perturbative evaluation of the generator functional of Eq. (7) implies a strong limitation. The naive perturba- 
tion expansion relies on the translation invariance of the vacuum, allowing weak inhomogeneous components in the 
external sources. This is enough to support sufficiently inhomogeneous field expectation value profiles when there is 
no gap in the excitation spectrum above the vacuum, such as in the case for a system of charges at finite density. 
But once the localized states arise through a gap in the spectrum then the weak external inhomogeneities are not 
sufficient to produce localized states on the mass-shell. In fact, all we achieve at finite temperature and vanishing 
charge density is a weak inhomogeneous vacuum polarization by virtual pairs. The creation of states with few on-shell 
localized charges remains a nonperturbative issue in this formalism [20]. 



III. FLUCTUATIONS AND DRESSING 



Both the external, classical source j and the dynamical current (0|j|0) = J generate vacuum polarization. The 
simple structure of the QED action allows us to compare these two polarizations. To this end we consider the 
expectation value of the EM field induced by weak external and dynamical, quantum currents. We shall seek the 
linearized expression for the EM field. The result will be equally valid for the transition amplitude formalism or for 
any double-time functional. The integration over the fermion fields in the generator functional (8) yields 



e iW[a,j] _ j D ^ e iW c [a-eaA] + ±AD 1 A+ij A 



where a denotes flipping the sign of the — component of CTP doublets 



A+ 
-A~ 



the generator functional of non-interacting fermions is denoted by 

W c [a]=Y,W n [a n ] 



(26) 



(27) 



(28) 



with 



W n [a] = -zTrln 



A-i, [>+ 

0n + I jr 



(29) 



being the generator functional for the connected Green functions of the current for the flavor n in the absence 
of electromagnetic interactions. Note that there are connected Green functions of arbitrarily high order for the 
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noninteracting Dirac-see because the current is a composite operator. We perform the shift of the integral variable, 
A — > A — D j, in other words, we use the equation of motion for the EM field to write 



W[&3] _ e -fj£>o~j J D ^ e iW c [a-ea(A-D j)] + iAD Q M 
— p iW[a+eaD j]-^jD j 



(30) 



with 



W[a] = W[a,j = 0] (31) 
being the generator functional for the current in full QED. Therefore, we obtain the relation 

W[a, j] =W{a + eaD j] - \jD j, (32) 

the reduction of the two- variable generator functional into a single variable functional. The simple form of the initial 
Lagrangian allows us to separate the tree-level dependence on the source j and to place the loop-induced dependence 
into the dressed fermion one-loop graphs. In doing so, we rely on the fact that the only part of the QED action which 
contains higher-than-quadratic terms in the fields is the minimal coupling. 

The natural use of such a reduction is the simplification of the equations of motion for the currents 

Jn = 5W ^t e ^ ] (33) 

Sa n 

and the EM field 

i a „ SW\a + eaD j] ~ ~ , . 

A = V e n D & l — °J± - D j (34) 

which can be written as 

A = Do(ae tr J — j) (35) 

where e denotes a flavor column vector made up of the electric charges, e tr = (eo, ei, • • • ). 

The external sources a n and j are the independent variables in these equations. When the functional Legendre 
transformation is performed then J n and A become independent variables but the equations remain valid by consid- 
ering j and a n as dependent variables. Comparing Eq. (35) with the second equation in Eqs. (18) the exact result 
[24] 

T[J, i] - T mech [j] + ^-Ab^A - AaJ^enJn (36) 

n 

follows. Here the first term represents the mechanical contribution of the Dirac-see and plays the role of an integration 
"constant" . There are two sources of complexities in the dynamics, the implications of the Pauli exclusion principle for 
the electric current as a composite operator and the electromagnetic interactions. The latter has a small parameter, ft. 
or e 2 to organize a systematical approximation scheme but the former type of dynamical correlations admit no small 
parameter and induce arbitrary high order connected Green functions for the non-interacting Dirac-see in Eq. (29) 
and arbitrary high order vertices in the effective action (36). 

The lesson of Eq. (35) is rather surprising, it is the absence of renormalization in the nonmechanical part of the 
theory. We shall see below, when a similar result is recovered in the loop-expansion, that this triviality results from 
a cancellation between the two channels the source j can induce EM field. What is observed here is that the relation 
among the expectation values and the external current j is not renormalized because the former already includes all 
the dressing. The nontrivial source of dressing according to Eq. (32) arises from the Green functions for the currents, 
the structure of the functional W[a] of Eq. (31). The j dependence comes in a trivial manner, dictated by the minimal 
coupling, gauge invariance. 

It is instructive to follow what happens with the vacuum polarization when we select the flavor n — as our valence 
charge, e v = eo, and cancel the external source a„ b n& ^ for the remaining flavors, called background charges. Let 
us introduce the notation a v = clq, 

1 _ 6W[a + e&D j) 

■Jv — j (3 1 ) 

oa v 
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and eb will denote the background charge vector, ejj r = (0, ei,---). The polarization effects now arise from the 
dependence of the electric current of any flavor on the external sources a v and j and one can find these contributions 
by expanding Eq. (33) in a v and in j. We have to keep in mind that the current of the valence charge, given by Eq. 
(37), includes all external source dependence and it is an independent variables of the effective action therefore, its 
dependence on j has already been accounted for. Thus the expansion can be restricted to the background charges 
which yields, after inserting it into Eq. (34), 



D 



c C V J V 



3 + vef 



S 2 W[a] 
SaSa v | S=0 



8 2 W[a] 



SaSa |a=o 



e b aD j 



+ 



(38) 



up to quadratic terms in the external sources for charge conjugation invariant vacuum. The first two terms on the 
right-hand side stand for the direct, tree-level part of Eq. (35) and the polarization effects are represented by the third 
and fourth terms in the linearized equation of motion. Recall that the external field a v is used only to generate the 
valence current J v . Therefore, a v in the third term should be expressed in terms of J v when the relation between the 
valence current and the EM field is sought. Thus this term represents the polarization effects due to the dependence 
of the background currents on the valence current. The last term shows the polarization due to the dependence of the 
background currents on the external current j. These two terms are not directly incorporated in the valence current 
dynamics and they represent the dressing in the equation of motion when a v is expressed in terms of J v . In other 
words, the dressing in the Maxwell equation (38) is due to the fact that our experimentally monitored valence current 
does not cover all charges in the system. It is completely natural that the uncontrolled charges generate the dressing, 
which gives rise to the proportionality between the averaged and the local EM field. 

IV. EFFECTIVE ACTION 

After some general remarks we now turn to explicit loop-expansion expressions for the quadratic part of the 
functionals mentioned in the previous Section. 

A. Connected Green functions 

The generator functionals of Eq. (26) for a single valence charge n = can be written after carrying out the 
integration over the fermion fields as 



JW[a,j 



Sj] = f D ^ e T. n W n [So, n d]+ijA+iAD- 1 A+0(a 2 ) + 0(A 3 ) ^ 

where W n is defined by Eq. (29) and a homogeneous, classical background charge is assumed to cancel the O(A) 
tadpole terms when the vacuum is charged. The improved photon propagator 



^=(^ - 1 -fl)- 1 
contains the one-loop polarization from all charges, 

ft = ae tr Gecr, 

where 

The background charge improved photon propagator, used later is given by 



D h = 



1 



(40) 

(41) 
(42) 

(43) 



including the polarization 



n 6 = ael r Ge b a. 



(44) 
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The CTP expression for the functional M^[a,j] is obtained when the photon field trajectories are closed in the 
functional integral of Eq. (39). We shall need the OTP result which is obtained by carrying out this functional 
integral for uncorrelated A ± trajectories, having different end points. The functional integration yields 



WM] - Wk T 




+ <D{h 2 ) +0 (source 3 ) 



(45) 



where G v = G 00 , G e = (Ge) . 

The dressed propagators are the sum of the products DxaD^a • • • oD n where Dj has the form shown in Eq. (12). 
The product displays the same structure as the factors Dj, namely the retarded or advanced part of the product is 
the product of the retarded or advanced parts, 



(£> lt xD 2t x • • • aD n a) a = D\D% ■ ■ ■ D a n , 



yielding 



(46) 



Dq _1 - n s 



(47) 



We need the real part of W, considered in real space as opposed to momentum space, for the effective action in 
Eq. (19). To simplify matters we exclude pair creation processes by restricting ourselves external sources with modes 
co, k < m and the heat and particle baths are chosen to be nonrelativistic, T, kp <C m. The one-loop expressions for 

both JRG and SG are vanishing on the photon mass-shell, the support of 3Z?o, according to the calculation reported 
in Section V. These properties allows us to simplify the real part of the products in the matrix elements in Eq. (45) 
and one finds 



3W[a,j] 



;(a,j) 



3?G„ e v $t(G v )a$t(D ) 

K e v M(b )a$t(& v ) $tb + R(D )R(f[)R(Do) ) 



(48) 



The form (12) of G and D gives after a lengthy but straightforward calculation 



KW[a,j] 





/ -kG" 


Ke v G n D% 


-G r 


e v G r D r \ 


I (a,j,a,j) 


Ke v DJyG n 


— K,D n 


e v D r G r 


_ D r 


2 


-G a 


e v G a D a 










K e v D%G a 


_ D a 





o / 



(49) 



where G a denotes (iRG) a . This result leads to the linearized expressions 

J = KG n a-Ke v G n D%]+G r a-e v G r Dlj 

A = - K e v D a G n a + KD n j-e v D r G r a + D r j 

J a = G a a-e v G a D%j 

A a = -e v D%G a a + D a j 



(50) 



for the expectation values in terms of the external sources for all charges or for a single valence charge, respectively. 

Notice that the physical expectation values J and A contain the retarded field of the physical sources a and j. This 
is not a dynamical breakdown of the time reversal invariance, but rather a trivial result of the boundary conditions. 
In fact, the open ended, null boundary condition at the final time leads to a destructive interference between the two 
time axes which cancels the advanced field generated by the physical sources [20] . The book-keeping variables a and j 
generate an O(k) time reversal invariant near field, as a result of the equal coupling of these sources to the dynamical 
variables of the two time axes by the K-dependent term in Eq. (14). 

After the formal manipulations we set a = j = to regain the physical case with unitary time evolution. The 
auxiliary fields are vanishing in this case, cf. Eq. (13), and the Maxwell equation, the expression of the EM field in 
terms of the dressed, retarded analogy of the Lienard-Wiecher potential, reads 



A = D r b j - e v D r J 



(51) 



where D^ 1 = D^~ l - nf . 
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TABLE I: Properties of the near and far propagator on and off the mass-shell 



Operator Mass-shell Off mass-shell 

UD n 1 

D% -D" 1 

Do ^0 

d^dI 1 - 1 1 

dI~ 1 D%D%~ 1 D^- 1 



B. Legendre transform 



The variation principle governing the expectation values is based on the effective action obtained by a functional 
Legendre transformation of SRW[a,j]. We mention first a peculiar feature of the loop-expansion. The independent 
variables of the generator functional for the connected Green functions W[a, j] are the external sources a and j which 
are classical, 0(h ) quantities. The result is the well-known equivalence of the expansion in h and in the number of 
loops in the Feynman graphs in the connected Green functions. But in the present case the electric current J is 0(K) 
and the expansion in h mixes different loop-orders in the effective action. 

The Legendre transform of a quadratic functional remains quadratic, and the two kernels are the inverse of each 
other up to a sign. The inverse of the block matrices of the quadratic functionals of Eq. (45) can be obtained by 
performing the change of variable of the Legendre transformation explicitly with the result 



--(j,A,.r,A a ) 



( 





Qa-1 








Cy 




Cy 


-KG r - 1 G n G a - 1 


V e v 


Dl- 1 


-Ke^PonG^- 1 + P of f) 



Cy 

or 1 

-K(P off + G r - 1 G n P on )e v 




(52) 



The operator P on projects on to the photon mass-shell and P ff = 1 — P on - Some important properties of D n 
and D' related to the mass-shell and used in the derivation are listed in Table I. For instance, the replacements 



0, D r 



-1 rjn n«- 



D^Poff and n^l-DJDg" 1 ) 



n^P on have been carried out 



in obtaining the 0{A a ) term in the effective action. The corresponding linearized, one-loop equations of motion are 



= e v A a + G a -\j a 
= Df- 1 A a + e v J a 



G^J-kG 7 *- 1 
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di- x a - K (n^p ( 



QUQa-lja + eA _ K (f> o// + G^G" P on )e v A a 

D*- l P off )A a - Ke v (P on G n G a - x + P off )J a + e v J 



(53) 



Notice that the auxiliary fields couple in a different manner to the physical external sources on and off the mass shell. 
The last expression is the vacuum polarization corrected Maxwell equation. It corresponds to the equation of motion 
Eq. (4.18) in Rcf. [21] when the effects of the physical sources which drive the system adiabatically to the desired 
initial state are incorporated implicitly in the propagators and the self-energy. 

The auxiliary fields are not needed for the dynamics of the observable expectation values therefore, it is natural 
to seek the effective action involving the physical fields alone. To this end we consider a and j as parameters in the 
generator functional of Eq. (49) and perform the functional Legendre transformation on a and j only, with the result 



kT[J, A] 




+ Ja'l + Aj", 



(54) 



where the inverse of the near field propagators is defined to be zero on the mass-shell. One may even keep only one 
field performing the Legendre transformation for only a single source with the result 



K T[A} = ±ADr 1 A + Af, 



(55) 



and 



kT[J] = -JG n ~ 1 J+Ja'. 



(56) 
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The source terms a" ,j" ', a' and j' are the sum of the action of different retarded propagators on a and j, their detailed 
form is such that Eqs. (50) are satisfied. Notice the need for n ^ to arrive at any variational scheme for the physical 
variables only [20]. 

What is the scale regime where the renormalization group treatment, performed in its complexity is supposed 
to preserve the form of these equations? The UV cutoff of the equations is the Compton wavelength because the 
vacuum polarization effects suppress the variation in the space-time within this distance scale. The infrared limitation 
is provided by the typical length scale of the state generated by the external sources, such as the de Broglie wave 
length. In fact, the smearing of the configuration beyond this length scale modifies the space-time dependence and the 
equations of motion. Once the blocking scale passes the quantum-classical crossover one arrives at the true, classical 
equations supported by decoherence. 



C. Radiation field 



The external sources or boundary conditions in time determine the radiation field in classical electrodynamics. The 
non-trivial issue here is the realization of this well-understood circumstance within the variational principle. The point 
is that the radiation field enters, in a singular manner, into the variational scheme of the actions (54)-(55). The first 
sign of this complication is the vanishing of the near field propagator on the photon mass-shell due to the principal 
value prescription in its Fourier integral representation. Another way to see that the radiation field does not enter 
into the quadratic part of the action is to recall the remark made after Eq. (12), namely that the far field propagator 
is antisymmetric thus the radiation field drops out entirely from the 0(A 2 ) part of the action. This state of affairs 
is in agreement with the formal time reversal symmetry of the elementary processes because the field operators act 
forward and backward in time by annihilating and creating elementary excitations. 

The retarded propagators may appear in the O(Aj') terms of the action and are attached to the external sources. 
The 0(A 2 ) part of the action controls the restoring force acting on the oscillation around the stable vacuum thus the 
external source induces a singular response on the mass shell, rendering the mass-shell radiation field nondynamical, 
to be settled by the boundary conditions in time. In fact, a simple way to arrive at the desired equations of motion 
without separating explicitly the on- and off-shell components of the equation as in Eqs. (53) is to write the source 
in Eq. (55) as j' = D n ~ 1 A on and note that the equation of motion A = A on contains the mass-shell component with 
a formal 0/0 coefficient. This scheme is actually realized by the usual m 2 — > to 2 — ie prescription for the free action 
when applied in classical field theory. 

When all fields arising from the reduplication of the degrees of freedom inherent in the double-time formalism are 
retained then the physical fields can be coupled to nonsymmetrical kernels. The result is the simple and natural 
preservation of the retarded and advanced solutions in the variational scheme. Such a treatment of the boundary 
conditions in time opens up possibilities to address the genuine, dynamical breakdown of time reversal invariance. For 
instance, the friction forces signal in classical physics that the coupling of our system to its environment breaks the 
time reversal invariance in a dynamical manner. Such kind of forces can easily be encoded in the variation principle 
of the double-time formalism. Another example is the dynamical building up of decoherence in a system due to a 
coupling to its environment with a gapless spectrum, to be seen explicitly below. 



V. ONE-LOOP POLARIZATION TENSOR 



We return in this Section to the loop expansion and present the photon-polarization tensor in the one-loop approx- 
imation. We shall work with the Fourier transform 



Jq 

where the integration in the Fourier space is denoted by 

The real and imaginary parts of the two-point function are defined in the space-time. Their Fourier transform is 

(itd) g = ±(d g + d*_ g ), (3(3% = 1(4 - &*_ q ). (59) 
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Our goal is to obtain 



ih 



(60) 



for a single charge in terms of the CTP propagator (B18). Because of the structure (12) it is sufficient to find G ++ 
and G + ~ . It will be useful to recall that this satisfies the Ward identity 



fc"GW = 0. 



(61) 



For the sake of simplicity we seek the detailed expressions below the pair creation threshold, q 2 < 4m 2 , and keep the 
density and the temperature low enough to keep the charges nonrelativistic. 

The CTP propagators are the sum of the contributions from the vacuum and from the environment. The latter is 
linear in the Fermi-Dirac distribution function n p — l/{exp(/3[e p — siga{p + 1}. As a result we have the similar 
separation in the polarization tensor (60), which will be the sum of three kinds of terms. The vacuum contributions 
are independent of n. The 0(n) terms contain a particle from the environment on the mass-shell. Finally, there will 
be 0(n 2 ) terms where both particles come from the environment and are on the mass-shell. 

As far as the environment is concerned, we shall have either finite temperature at vanishing density or finite 

density at vanishing temperature. We write the two-point function as G — G vac + G env where G env depends on the 
temperature T or the Fermi momentum hp. The matrix element H — h of the polarization tensor in the vacuum is well 
known, its form is 



q++i»> = "2 




1 4m 2 / 4m 2 
— - 1 arccoti / — 1 — 1 



(62) 



with a = e 2 /47r. The matrix element 



w vac,q 



-128n 3 iam 2 ( g» v 



2m/ 



f 5{q 2 + 2pq)S(p 2 - m 2 )e(-p° - q°)e(p°) (63) 
Jp 



is vanishing below the pair creation threshold because it consists of on-shcll amplitudes only. 

The environment breaks the formal Lorentz covariance of G env . Let us denote the temporal unit vector by 
which assumes the form = (1,0) in the rest frame of the environment where the expressions for the distribution 
functions of Eqs. (B19) are valid. We use two invariants q 2 = — \q — u(uq)] 2 and v = uq/\q\ = oj/\q\ with q^ = (u>, q) to 
parameterize the momentum dependence. For each CTP index, covariance with respect to three-dimensional rotations 
suggests that we use three three-dimensional scalars whose number is further reduced to two by gauge invariance. 

These two sets of scalar are the easiest to obtain in terms of A = g^G"^ and B = uGu. The final useful form is 



(G g r=B q 




- 2 [B q (l 



u 2 )-A q ] 




,0 T, 



(64) 



where the space-time tensor structure is explicitly shown for the environment contributions with n = q/\q\, L = n®n, 
T=l-L. 

The expression (64) makes the space-time tensor structure explicit. The CTP index structure is contained in the 
Lorentz-scalar CTP matrices A and B which will be detailed now. These matrices are the sum of the vacuum and 
environment contributions, 



Aq X vac _q + A, 



env,q 



(65) 



where A stands for A or B. As far as the vacuum contribution is concerned in the framework of the gradient expansion 
it contains G++^ only which gives 



a T „ 2 f 1 ) _ (SRAac)q 



(ftB vac ) q - — ^q 2 



-1 3(1 -v 2 ) 



(66) 
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and QAvac = ^B vac = with the help of Eq. (62). The structure of Eq. (12) is preserved by the environment 
contribution, and we find the real and imaginary parts 

= \ f 27rS (l 2 + 2pq)2ir5(p 2 - m 2 )sign(/ + q°)(n p + n_ p )X p , q 

(QX env ) q = i J 2i:5{q 2 + 2pq)2ir5(p 2 - m 2 )[n p (n q+p - 1) + n_ p (n_ p _ 9 - l)]X Ptq , (67) 

where P denotes the principal value integral and the kinematical factors are 

A Pt g = 8(m 2 -pq), 

B p , q = 8(p 02 +p q -±p q y (68) 

We shall consider the cases of finite temperature and vanishing density and finite density and vanishing temperature 
when the one-particle distribution function n p generates the scales k c h — T and k c h — kp, respectively and keep the 
environment in the nonrelativistic regime, T,kp <C m for simplicity. The frequency integrals are carried out in Eqs. 
(67) by means of the residuum theorem, resulting in three-dimensional integration over the momentum p. By splitting 
these integrals into two parts, corresponding to \p\ < m and \p\ > m and by finding an upper bound for the latter 
contributions one can verify that the simple nonrelativistic expressions give reliable leading order approximations. 

What are the kinematical regions contributing to the various pieces of the two-point function for a nonrelativistic 
environment? At low temperature and vanishing charge density it is enough to retain the 0(n) terms which decreases 
rapidly when \q\ is increased. The situation changes drastically at finite density because the gap disappears in the 
excitation spectrum. 3?X ++ , given in the first equation of Eq. (67) describes a particle- hole pair where one member 
of the pair belongs to the environment and is on the mass shell. It is advantageous to use the dimensionless ratio 

r - q \^: Pq ° (69) 
2\q\k ch K ' 



which gives 



(g + rk ch n) 2 (rk ch ) 2 

w ~ ^ ^ (70) 

2m 2m v ' 



indicating that the Pauli blocking allows particle- hole excitations at small q for |r| w 1. Thus $tX ++ is expected to 
take larger values for q <C m at \r\ rts 1 only. Both the particle and the hole are on the mass-shell in the remaining 
expressions in Eqs. (67) because the time ordering is the only mechanism to produce off-shell amplitudes in the 
CTP propagators. The distribution function n p restricts the integration approximately to \p\ < k cr . The conditions 
p 2 = (p — q) 2 = m 2 for nonrelativistic four momentum p with p° w m give rk cr = pn, the integration is over a plan, 
orthogonal to q. The largest area corresponds to r = and the remaining functions in Eqs. (67) tend to be large for 
r w 0. 

The actual calculation of the integrals leads to the expressions 

= X+ + X+ „ , 

Q q u ,q — q u ,q' 

Q V — a". a o u .O'" 



9 v -q",q q",q> 

= X % q0q + X z _ q0q , (71) 
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where the case of finite temperature and vanishing density gives 



A + = 



Aam 2 T 



A- = 



A 1 = - 



Bt = - 

1 



BZ 



Bl 



Aam 2 T 

\q\ 

Aam 2 T 
Aam 2 T 

7r|q| 

Aam 2 T 
Aam 2 T 



2m 2 J 

2 

_ 



AaT 2 ~ 



2m 2 
2m 2 

„2 



K, 



1 + 



Am z m 



2aT 2 ~ 



Am 2 

„2 



— J, 



1 + 



Q 



q_ 

m 



, 2 +- 

Am z m 



K, 



(72) 



where 



f°° dzz 
Jo u z 



In 



r z + z 



r z - z 



Jo = 



^-n z Q(l - r z )sign(sign(g°)o) 2 

LU Z 



dzz 



n,0(l-r) 



(73) 



U z l + e -V^ 2 +2 2 + (/39 ) 2 +2/3 9 ^ ' 

with r z = 2 ^| T + <1> Z = ^/ ((3m) 2 + z 2 and n z = l/( e V(/ 3m ) 2 + z2 -f 1). The results for finite density and vanishing 



temperature are 



A + 



Bt 

q 



2ak F m 
ir\q\ 

ak 2 F m 



2akpm 
n\q\ 



b: 



A 1 

■*q 



\Q\ 

k 2 F i 



ak 2 F m 



2m 2 

„2 



2m 2 J 

2 



M, 



9' 



1 

1 + 

1 + 

1 + 



Q 



Am 2 
Am 2 
2m 2 J 

r2 



Am 2 



(£_ 

rn 
<t_ 
m 

N q , 
m 



M a 



N„ 



(74) 



with 



L q = r+i(l-r 2 )ln 



1 



r - 1 

M q = e(l-\r\)(l-r 2 ), 

'l-r 2 \q\>2k F , 
l-r 2 \q\<2k F , 

k -i |g|<2fc F , 



N a = < 



- 1 < r < 1 



LgJ 

2k f 



Mi 

fe F 



< r < 1 



< r < 1 



Mi 

fc F 



(75) 



The variable r in these equations is given by Eq. (69) after the replacement k c h = k F and uj p = m. 

L is the relativistic generalization of the Lindhardt function. The last two expressions in Eqs. (74) give the 
relativistic version of the decoherence factor of the Coulomb field, considered in Ref. [22]. The expressions (72)-(75) 
are valid for nonrclativistic electron gas but naturally comprise the relativistic kinematics for photons. 
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VI. REDUCED DENSITY MATRIX 

We use now the functional in Eq. (8) to obtain the reduced density matrix for the EM field. A single charged field 
is considered for simplicity and it will be treated by CTP boundary conditions and the OTP boundary conditions are 
used for the EM field. The generator functional in Eq. (8) with a = j = 0, considered as a functional of the final field 
configurations gives the reduced density matrix 

p[Aj x ,AjJ = J D[i](n^[^]^[4])e i ^^^on 1 -en^„ + iAA - 1 A_ (?6) 
The integration over the charged field can be carried out with the result 




where the effective bare action is 

S eff [A] = l -A{D^ U)A + 0(i 3 ), (78) 

II being given by Eq. (41). This functional integral contains information about the dynamics described by the density 
matrix in a manner similar the path integral describes the transition amplitude between pure states in the single time 
formalism. The real part 3f5L4] is identical to T[A], obtained as the Legendre transform of ^R.W[a, j] for a = 0, it 
determines the expectation value of the EM field and will be used to extract the polarization induced by the charges. 
The imaginary part controls the width of the peak in the reduced density matrix, and the decoherence and consistency 
of the EM field. 

The normal modes of the quadratic part of the action in the path integral expression for the transition amplitudes 
are plane waves. It is easy to see that H q is not diagonalizable in the time axis index only. The normal modes of the 
CTP action (78) couple the wave vector q = (oj, q) with the CTP index and are labeled by the wave vector, space-time 
and CTP indices. The dynamical role of each normal mode is characterized by three numbers, 

n?"" = -/ dxe^mpwx]), 

11^ = - J dze^RiOm 

nr - / dxe^iiTU^}), (79) 

by applying the parametrization (12) for the current-current two-point function G. It is easy to check that the Fourier 
transforms II" and II 1 are real and n/ is purely imaginary. In the one-loop approximation the real quantities arise 
from the interactions within a single time axis and the purely imaginary quantity corresponds to correlation between 
the time axes. This correlation, realized by Green functions connecting the two time axes, arises because the charged 
particles which are exchanged in the process represented by the Green functions are entangled with the EM field. 

We shall consider these quantities in the (q, ui) space. The zeros of II™ locate the mass-shell of the quasiparticles, 

represented by the normal modes of the diagonal, H — h or blocks of the quadratic bare action. The imaginary 

part of the two-point function in the space-time, which is given by ill^ in Fourier-space, controls two, superficially 
different dynamical processes. On the one hand, the decay of the quasiparticles, the inverse life-time of a quasiparticle 
can be identified by 11^ evaluated on the mass-shell. On the other hand, according to the structure displayed in Eq. 
(12), this parameterizes the suppression of the off-diagonal matrix elements of the reduced density matrix, realized by 
the mode ui, q of the path integral. Thus the consistency and decoherence of the quasiparticle modes on the mass-shell 
have the same dynamical origin as the finite life-time of the quasiparticles. This is the expected relation between the 
dynamical breakdown of the time reversal invariance and the classical limit. 

We now turn to the discussion of the qualitative features of 11™ and 11^ in describing the polarizability of the charges 
and the impact of IP on the classical limit. 

A. Electric and magnetic susceptibilities 

To make contact with the usual three-dimensional notation we introduce the parametrization A^ — ((f), A) of the 
vector potential, giving rise to E — — V</> — 8qA and B = V x A. We consider quadratic action in the EM field in a 
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formal manner, as a book-keeping device for the equation of motion (51), namely we assume that the actions contain 
the inverse of the retarded Green function. On the one hand, we can introduce the retarded electric and magnetic 
susceptibilities \ an d X, respectively with the phenomenological parametrization 



MA] = \ I ( 

Jx, y 



Ex^XyEy B X jl X yBy 

) (80) 



of the Maxwell action. On the other hand, the effective action for the EM field expectation value yielding the equation 
of motion is 

T = ^AD r ~ l A (81) 

where D r ~ x = D^ 1 — e(G ++ — G H )e according to Eq. (47). Therefore the equivalence of the two quadratic form, 
expressed by the equation 






1 -nv \ ,_ . 
-uv vH-T )-^- e) T, 



(82) 

follows in Fourier space with X r = X ++ — X + . The left hand side contains the free inverse propagator and the 
retarded self energy as given by Eqs. (64) and (66)-(67). The right-hand side is S 2 So[A]/5A- q SA q with v — w/|q|, 
n = q/\q\, L = n®n, and T = 1 — L. We insert into the right-hand side the expressions e = 1 + \ an d ft = 1 + x 
which gives, after canceling the free inverse propagator on the left-hand side by the contributions of the one on the 
right-hand side 



M ~ q 2 ' 

x q - gl 2q J " (83) 

Let us now see the more detailed, numerical features of the polarization functions at finite density, obtained by 
using Eqs. (71), (74)-(75) to compute X r q = X+ + — X+~ in Eqs. (83). The susceptibilities display a divergent 
structure which is well known for nonrelativistic Coulomb gas. The magnetic sector and the relativistic effects for the 
EM field are retained in this calculation which slightly modifies the picture. The susceptibilities, q 2 e u ,q and q 2 ftu:,q 
are shown in Figs. 1 for m = 1 and fi = 0.1. They all take appreciable values for |r| w 1, cf. Eqs. (74)-(75) and 
(70). The quasiparticle conditions Sfie^q = 0, $tft. u ,q = have solution along the valley of the line r = 1, shown 
on the (|q|,w) plane in Fig. 2. The approximately linear parts in the electric and magnetic sectors belong to the 
zero sound, the back turning sections correspond to the plasmonlike excitations. The EM field may have relativistic 
energy and momenta and the relativistic correction factors 0(q 2 /2m 2 ) in Eqs. (74) modify the plasmon lines in an 
essential manner. They prevent the plasmons from becoming long range and push the frequency up sharply as the 
wave vector tends to zero in the electric sector, as seen in Fig. 2 (a). For w > 2 we naturally run into the pair 
creation singularities. The normal modes of the magnetic sector behave in a rather peculiar manner, the plasmon 
line, shown in Fig. 2 (c) follows approximately a nonrelativistic dispersion relation curve of mass m e // = 10 -3 even 
in the relativistic domain. 

The screening of the electric sector manifests itself in the nonvanishing value of the product q 2 $le in the limit q — > 
and the resulting infrared divergence in |e| as shown in Fig. 3 (a). The valley in |e| identifies the strongly coupled 
modes, they are along the plasmon line of Fig. 2 (a) and on the continuation of the zero-sound line. The zero-sound 
line shown in Fig. 2 (a) is weakly coupled due to the short life-time. The magnetic sector is qualitatively similar, \ft\ 
plotted in Fig. 3 (b) indicates strongly coupled magnetic plasmons and weakly coupled magnetic zero sound modes. 
The magnetic field is screened as well as shown by the nonvanishing value of q^Rft in Fig. 1 (c) in the infrared, q — > 
limit. 

The physics is simpler at low temperature and vanishing density due to the gap in the excitation spectrum for 
the charges. One finds qualitatively similar frequency and wave vector dependence as displayed in the figures above 
except that the overall scale is suppressed by exp(— m/T). 



B. Consistency and decoherence 



The reduced density matrix, given by Eq. (76) can always be diagonalized and the resulting basis consists of 
completely decohered states. But one can gain no more insight into the classical limit by using this basis. Instead, we 
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FIG. 1: The susceptibilities, (a):g 2 3^x = <? 2 (9fe — 1), (b):q 2 ^e, (c):q 2 $lx = q 2 ffifi> — 1), (d):</ 2 £y/2 as functions of uj and \q\. 




0.01 0.02 0.03 0.04 



(a) 




0.001 0.002 0.003 0.004 0.005 0.006 



(b) 



0.01 0.02 0.03 0.04 0.05 



(c) 



FIG. 2: Quasi-particles on the (|qr|,w) plane, (a): In the electric sector they are the roots of the equations = 0. The zero 
sound and the plasmon modes correspond to the straight line starting at the origin and after the turning point, respectively. 
The magnetic quasiparticles are at 5R/i = 0, (b) is around the origin and (c) is on a larger part of the plane. 



are interested in the role of physically motivated states which persist at both side of the quantum-classical transition. 
These correspond to quasiparticles, denned by the 0(A 2 ) part of the bare action (78). Therefore, we seek the 
contribution of the quasiparticles to decoherence in what follows. It is easy to see by means of the structure (12) of 

G and the expression (41) for the self energy that the imaginary part of the effective bare action of Eq. (78) 
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depends on the combination A + — A of the EM field only. Eq. (82) can be used again to parameterize it as 

1 



SiS[A] = - I (E d e d E d - B d fi d B d ) (84) 



with E d = E + — E , B d = B + — B where the decoherence parameters 

ni 

d _ u g 

Q q 2 

. Bin - iv 2 ) - Ai 
K - -± — -r 1 q -, (85) 

are expressed in terms of the imaginary parts introduced in Eq. (71) and are plotted in Fig. 4. 535 [A] governs 
the suppression of each plane wave mode along the off-diagonal direction in the reduced density matrix. It is easy 
to see that this functional is semidefinitive because it measures the phase space available for particle-antiparticle or 
particle-hole excitations which incorporate decoherence. Its value, taken on the quasiparticle line of Fig. 2, yields 
the inverse life-time. We see that the zero-sound quasiparticles are rendered unimportant by their short life-time. 
Naturally, higher loop contributions to the effective action generate imaginary part and finite life-time beyond the 
zero-sound line by the multi particle-antiparticle or particle-hole excitations, but this remains a weak effect. 

When considered from the point of view of the quantum-classical transition, the effective bare action monitors the 
impact of the entanglement between the EM field and the charges. One measure of this entanglement is decoherence 
which corresponds to an instantaneous state. But we find more information in this action. The suppression of the 
contributions in the path integral representation of the density matrix as the functional of the difference of the two 
field trajectories carries information about the consistency of histories of the EM field [3] . The contribution of a plane 
wave, a normal mode of the one-loop dynamics of the EM field, to the conditional probability distributions is closer 
to the one expected in the standard, classical probability if 95 [A] increases with A + — A~ because the quantum 
interference term which violates the additivity of probabilities is suppressed by e _ss [23]. 

The lesson is that the classical probabilities are recovered by plane wave modes whose dispersion relation is that 
of the short life-time quasiparticles. Their short decay time assures the dynamical breakdown of the time reversal 
invariance, needed for the classical limit. These modes which maximize 95 are located around the long wavelength 
part of the zero sound curve at finite density and vanishing temperature. The suppression of the off-diagonal elements 
of the reduced density matrix for these plane wave EM fields shows that these field configurations become good pointer 
states in the classical limit. 



VII. CONCLUSIONS 



The problem of deriving classical electrodynamics from QED is pursued further in this work. It is based on the 
monitoring of subclassical fields, and the expectation values of local operators with space-time resolution limited 
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formally by the UV cutoff of QED. The final goal is to weaken the resolution and thereby recover the macroscopic, 
classical dynamics for the space-time dependence of the expectation values. As a preparatory step, we make a coarse 
graining by eliminating the charged degrees of freedom with a fixed, microscopic space-time resolution. Polarization 
effects, arising in this manner, reflect competing coherent and decoherent phenomena, both at the elementary and at 
the collective level. 

The double-time axes formalism used in this work was motivated by the implementation of initial conditions, as 
opposed to the boundary value problems within the variational method. This step leads to the double-time axis 
formulation in a natural manner. The double-time axes formalism is well suited to follow the eventual breakdown 
of time reversal invariance at the quantum level, driven dynamically by degeneracies or externally, by the boundary 
conditions in time. The breakdown of the time reversal invariance appears at the level of the expectation value of the 
EM field as the radiation field. A simple microscopic way of separating the EM field into near and far components 
emerges as well for they represent pure state self-interactions for the charges and generate entanglement between 
charges and the EM field, respectively. 

Some well-known issues appear in a new point of view, such as the origin of the dynamics for the free-field components 
in the variational equations. The free-field component belongs to the null space of the kernel of the free action and drops 
out from the linearized equation of motion. This problem is avoided in the double-time formalism, and Feynman's it 
prescription seems to be the analogous treatment in the single time formalisms of both quantum and classical field 
theory. 

Another question which can easily be addressed is the relation between the classical limit and the dynamical 
breakdown of the time reversal invariance. The block structure of the CTP propagators shows that a shorter lifetime 
of the quasiparticles means more consistent EM field trajectories and a narrower peak in the absolute magnitude of 
the reduced density matrix \p{x, y)\ around the diagonal, x — y. Another way to see this relation by means of the EM 
field expectation value is to note that the purely imaginary far field component of the CTP polarization tensor, ±l£ q , 
plays a double role in the dynamics. On the one hand, it characterizes the mixed states arising from the charge-EM 
field entanglement, and on the other hand, it leads to a phase delay in the propagation of the EM field and thereby 
breaks the time reversal invariance. 

These results open the way for further developments, and we close by mentioning some of them. The subclassical 
fields are expectation values but they are not necessarily classical quantities because they may display the microscopic 
structure in space-time where decoherence is not efficient. In order to identify initial conditions which lead to classical 
fields, we need to control either the decoherence itself or the characteristic scale of the field expectation value. The 
latter is determined by the two-point functions and its handling requires the extension of the present scheme to 
external sources coupled to the bi-local product of the elementary fields. Such an explicit treatment of the correlation 
functions would, in addition, allows us to test different collapse scenarios of the wave function. 

Another extension of the present study is the detailed analysis of the effective action for the EM field and the 
electric current to derive a variational scheme for magnetohydrodynamics. The more accurate determination of the 
susceptibilities from first principles may lead to new experimental devices to identify some signatures of the underlying 
microscopic quantum theory by precision measurements in classical electrodynamics. 
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Appendix A: CTP propagators in the vacuum 

The CTP vacuum propagators for scalar, vector, and spin-half particles are summarized in this Appendix. The 
relation between expectation values and the parametrization of the propagator is achieved by path integration, and 
the explicit expression, for the different expectation values are computed in the operator formalism. 

1. Neutral bosons 

We introduce the propagator for a boson field (f> with or without vector indices by means of the free generator 
functional, 

e iw\j] = J Dfyetf-D-^Mt 

= e~^- brj (Al) 

which gives, in the Heisenberg representation, 

= ^(O||^){^|T[^ a & ]|O> = (O|T[^ a fe ]|O> =z J D++, 



i 



= £>|&|n)<n|0 o |O) = (0|^ o |0) = iD+ 

OIJa OIJ^ 

. S 2 W[j] 



b ' 



6 ' 



J2(Q\<t>a\n)(n\cf> b \0) = (0\cp a <p b \0) = iD~+, (A2) 



written as 



A D D+~\ / D -D-+A = ((0\T[<PM0) (0\MaO\0)\ (m 

\D-+ D—) ab \D-+ -£>+ ) ab \ (O|0 a fc |O) (0|T[^ a ]|0)*y 

by means of the notation D ++ = D. The basic CTP identity, 

T[cb a cb b ] + f[cj) a cb b } = <\>M + 6 O) (A4) 

written as 

D-D* = D~+ - D-+* (A5) 

allows us to parameterize the propagator in terms of three real functions, the near and far field Green functions, D n 
and , respectively, and the common imaginary part, D- L = as in Eq. (12). The detailed form of the real and 
imaginary parts, 

KD aa , = -±sign(i -t')i{0\ [0 , 4> a >]\0) = D2a„ 

D la a' = ~{0\{<f> a ,4> a '}\0), 

KD-J = -\{mM\V)=D f aa „ (A6) 
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justifies the introduction of retarded and advanced propagators 

aa 

Te(±(t-t'))i(0\[<f> a ,<l> a ,}\0). (A7) 



^aa' — ^aa' ^ ^ aa' 



The expectation values of the propagator (A3) are easiest to compute in the operator formalism by means of the 
quantum field 

<j>{x) = f [a(k)e- lkx + a){k)e lk '% (A8) 
Jk 



where the notation 

d 3 k 



1=1 — 

Jk J (27T 



) 3 2co k 



(A9) 



ujk = \/m 2 + k 2 has been introduced. The nonvanishing canonical commutation relation for the creation and annihi- 
lation operators, 

Hp), at(p')] = (2^) 3 2 Wp «5(p - p'), (A10) 

gives, in a trivial manner, 

iD^ix.x')) = f e ± ife (x-x') (An) 

Jk 

for a scalar particle. The multiplication by the Heaviside function, 

^ r g-iw(t-t') 



<d{t-t')iD-+{x,x')) = — / / 



e — iu h (t — t')+ik(x-x') 



e 



-iui(t— t')+ik{x — x') 



1 f e - iw (*'-*) 



1 r p —iLO(t -t) r 

G(t' -t)iD+-(x,x')) = — / / e -*-n(* -*)+*(« -«) 

r p —iu(t—t')+ik(x—x') 

= -i / (A12) 



gives the usual causal propagator and, we finally find 

( D D+-\ I -2mS(k 2 -m 2 )e(-k^ 
\D-+D-) k \-2ni8(k 2 - m 2 )6(fc°) -p^ 

at vanishing temperature. 

The Fourier transforms of the real and imaginary parts of the causal propagator arc 

Dl = P- ' 



(A13) 



k k 2 — m 2 ' 

Afc = -TrS(k 2 -m 2 ). (A14) 

In computing the real and imaginary parts of D~ , one has to keep in mind that they are defined in Eqs. (A6) in the 
real space-time as opposed to Fourier space, and the expression 

D-+ = -m [ e- lk{x - x 'h{k 2 - m 2 )[9(fc°) - 9(-fc )] - m [ e- lk(x - x ">S(k 2 - m 2 )[9(fc°) + e(-fc )] (A15) 
Jk Jk 

yields 

D{ = -TTi6(k 2 - m 2 )e(k°) (A16) 
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and 

D l = 9 \ ■ »nv ( A17 ) 

fc k 2 - m 2 ± iesign(fc°) v y 

Note the remarkable feature of the near field propagator D" in the first expression of Eqs. (A14). It is the only 
part of the CTP propagator which is nonvanishing off the mass shell, the support of the other functions Di and is 
just the mass shell. Furthermore, diverges in momentum space as we approach the mass shell, k 2 — > m 2 ; therefore 

D n-l D f = Q) ( A1 g) 

cf. Eq. (Dll) in [20] for the regulated expression. 

2. Charged bosons 

In the case of a non-Hermitian boson field the free generator functional is 

+** - fm^^^ - (ad 

and the CTP propagator block-matrix elements, identified in a manner analogous to Eqs. (A2), are 

A D D+~\ = /<0|T[^t]| ) <0|^ o |0) \ , , 

\D- + D-J \ (O|0„^|o> (O|T[0^t]| O )*y 

and they can be written as in Eq. (12) due to the CTP identity 

T[<f> x 4,]+T*[<f> x 4,] = 4,4> X + 4> X <I>1 

D-D* 1 = D-+ + D+-. (A21) 

It is the charge conjugation invariance which guarantees that the detailed form of the functions Di, D n and D> is 
the same as in the case of neutral particles. The non-Hcrmitian quantum field 

<j,(x) = f [a(k)e~ lkx + b\k)e lkx ] (A22) 
Jk 

involves the creation and annihilation operators whose nonvanishing canonical commutators are 

Hp), a\p')\ = [b(p), 6t(p')] = (27r) 3 2ov (A23) 

These operators transform under the charge conjugation transformation C, C 2 = 1, as C' f a(p)C — b(p), C^b(p)C = 
a(p) and C^tp x C = ip%.. Now, the relation D + ~ = D h actually establishes the equivalence of vacuum expectation 
values computed in the particle and the antiparticle sectors, and it originates from the charge conjugation invariance 
of the commutation relations (A23) and of the vacuum. When the propagator is expressed in terms of D and D '", 
the complex conjugation like in Eq. (A5) reverses the order of ip and ip> , and one needs another complex conjugation, 
this time embedded into the charge conjugation operator, to place the dagger at the appropriate operator in Z? H , 

iD+- = (0|^ o |0) = (0|CWtc|0>* = iD -+ c *. (A24) 

In the presence of a noncharge conjugate invariant vacuum, D A ^ D h . 

3. Fermions 

The generator functional for the propagator, 

e iW[j,j] = J D ^ D ^ e 4.G- l 4+fj4+4-3 = er^-G-o^ ( A25 ) 
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leads to 

A G G+-Y P = ((0\T[^]\0) (0|(^|0) \ , . 

\G-+G-) xy V-<0|^|0) (0|T[( 7 V y )^(^ 7 0)«]|0)*y 

by following the strategy of Eqs. (A2). The CTP identity, 

T[V>$ + f [^] = -^V + (A27) 

now leads to the parametrization 



G G+-\ = (G n + iGi Gf-id 
G~+ G—j ~ [-Gf-iGi -G n + iG, 



(A28) 



The field operator is written as 

%l>{x) = f — Y[b a {k)u^ a \k)e- tkx + dl{k)v {a \k)e tkx ] (A29) 

where {J: — m)u(k) = m)v(k) = and the creation and annihilation operators are defined by their nonvanishing 
anti-commutation relations 

{b a (p),bl(p')} = {d a {p),d\{p')} = {2nf^-8 a ^{p-p'). (A30) 

The charge conjugation is represented by C~ 1 ipC — Uip* with U = i"f 2 and the identity (A27) allows us to write 

i(G + ~)t r y = (U(0\C^ x ^ y C\0)*U) af3 = {UiG- +c U)* x f 
iG-~ a0 = (j iG yxl °r^. (A31) 

The Fourier decomposition (A29) gives the explicit expression 

(££)-<*♦..>(_£,-£:)• 

where the scalar propagator in the right hand side is given by Eq. (A13), in a standard manner. 

Appendix B: CTP propagators at finite temperature and density 

We extend the propagators recorded in the previous appendix for the presence of heat and particle reservoirs. 

1. Neutral bosons 

It is advantageous to use a finite quantization box, where the commutation relation (A10) reads 

[a{p),a\p')]=2Vw p S p , p/ (Bl) 

and the quantum field becomes 

m = 2~kr k EHfc)^ + a\k)e ik % (B2) 
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The relation 



Tr[e-^4>(x)4>(x')} = 



{n} k 



—ik(x—x / ) 

2Vuj k 



-({n}|e-^|{n}) 



k{x — x' ) _j_ g— ik(x— x')} 



){{n}\e^ H aHk)a(k)\{n}) 



I g-ifc(x-x') JJ^ e -flne h + ^_( e <fc(x-x') + e -*fc(*-*'))({ n }| e -0» n ( fc )|{„}) 

•'fe fen r„i !. Wfc 



-ik(x-x') 



-/3£fc 



En e " 



where i?o is the free Hamiltonian gives 

i£T + (a;,a;') = 



_ Tr[e-^0(x)0(^)] 



with -D T 1 being taken from Eq. (All) and 



AD = [ e 
Jk 



iD~+ (x, x') + AD(x, x'), 



tk ( x - x ') ^5(k 2 - m 2 ) 

e /3e fc _ 1 ' 



Similar steps produce 



iD + -(x,x') 



Tv[e-^ H <j){x')(j>{x)] 

Tr[e-P H ] 
iD^Z (x,x') +iAD(x,x'). 



The lesson is that finite temperature modifies the imaginary part of the propagator only, where 

2nS(k 2 - m 2 ) 



AD, 



(B3) 



(B4) 



(B5) 



(B6) 



(B7) 



and AD R = AD A = 0. The independence of the retarded and advanced propagators from the temperature can be 
understood by recalling that these propagators correspond to fixed initial or final conditions. 



2. Charged bosons 

The non-Hermitian field 

k 

involving the operators which satisfy the nonvanishing commutation relations 

[fl(p),flV)] = [*(p)>iV)] = 2Vw p J p ,p, 

give rise to 

-P(H-nQ) l 



Tr[e " ''-"'()(■'• 



)4t{x')\ = f 
Jk 



e ik(x-x') e -Pn(e k -n) 
k n 



(B8) 



(B9) 



(BIO) 



EE^(WI^ ( ^" Q) (* t (fe)i(fe)e tt( "' ) + aHk)a(k)e- lk{x - x ">)\{n}} 
i — - x — j 2Vu!k 



{«} * 



26 



which in turn gives 

iD-+(x,x') 



Tv[e-P( H -vQ)<j>(x)tf(x')} 



X r [ e -/3(ff-pQ)] 

- rU T ^ x) + ]^e ^ + y_ e ^ 

■ /^ /* -ifcfx-x'l /2^(fc 2 -m 2 )9(-fc°) 2^(fc 2 -m 2 )e(fc°)\ 

= *o + / fe c •*<» > ( + J • (B11) 

The modification due to the environment is the same in Tr[e~P( H ~ ,i Q' > $ (x')(j>(x)] as in Eq. (BIO), leading to the 
same conclusion as in the case of a neutral boson, namely that the only modification of the environment is D ±T = 
D^Zo + AD > and 

^(fc 2 -m 2 )6(-fc°) ^7n5(fc 2 -m 2 )e(fc ) 
lAAfe = ~ % e P(c h +,) - 1 1 e^(^-M) - 1 ' (B12) 

3. Fermions 

The field operator 

^(x) = 1 V — [b a {k)u^\k)e-' kx + 4(fc)u( a )(fc)e ifcx ] (B13) 

ka 

and the anti-commutation relations 

{b a (p),bl(p')} = {d a (p),dl(p')} = V^5 a ,p5 p ^ (B14) 



give 



{n} 



(d t (fc)d(fc>< Q Hfc)u (Q) (fcOe i * I_ifcV +6(fc)6 t (fc')« (a) (fc)« (Q) (fc')e ife ' a; ' _J ' £a; )|{«}) 
EE: 



, 2F 2 w 2 

{n} fe K 



({fi}|e-^-"°>((i t (fc)(i(fc)(it-m)e tt ( 1 - 1 '' + (fc)0& + m) e - ife ( x - x '))|{n}) 

= E E ^7 (* + ^K^-^ ({n}|e"^ ff -^ |{n}) £>0 
{«} fe 



+EE- 



m 



({n}\e-^ H -^(dHk)d(k)(fc-m)e lk(x - x '^ - b*(k)b(k)(]k+ m)e- lk{x - x '^)\{n}) 

= f{fi+m)e- tk{x - x "^\\{l + e- l) ^-^) 
Jk k 

k k 

This result, together with the analogous one for Tr[e _/3 (- H ~ M< 2)'i/5(x')'0( a; )] yield 

G ±T = -{ip+ m)(D±T Q _ AD), (B16) 
where the last term is given by Eq. (B12). 
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4. Summary 

The different propagators encountered in the discussion presented above can be summarized as 



D k = 
and 



k 2 — m' 2 -\~ie 



-2TTi5(k 2 - m 2 )6(-fc°) 



2m5{k 2 - m 2 )0(k°) 



i2TrS(k 2 - m 2 )[9(fc°)n+ + Q(-k°)nr] 



k 2 — m 2 —ie 



1 l 
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(B17) 



k' 2 —m' 2 -\-ie 



2TriS{k 2 - m 2 )e(-k Q ) 



2mS(k 2 - m 2 )6(fc°) 
where the occupation numbers 



k 2 —m 2 - 



nt = 



■ 2m5{k 2 - m 2 )[Q(k°)n+ + Q(-k°) 
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1 -1 

-1 1 , 
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e /3(e fc =FM) _ 1 



/ ~ e /3(e fc =FAi) _|_ 1 



(B19) 



are used. 
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